Self- Assembly of Glycine on Cu (001): the Tales of Polarity and Temperature 
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Glycine on Cu(OOl) is used as an example to illustrate the critical role of molecular polarity 
and finite temperature effect in self-assembly of biomolecules at a metal surface. A unified picture 
for glycine self-assembly on Cu(OOl) is derived based on full polarity compensation considerations, 
implemented as a generic rule. Temperature plays a non-trivial role: the ground-state structure at 
K is absent at room temperature, where intermolecular hydrogen bonding overweighs competing 
molecule-substrate interactions. The unique p{2x4) structure from the rule is proved as the most 
stable one by ab initio molecular dynamics at room temperature, and its STM images and anisotropic 
free-electron- like dispersion are in excellent agreement with experiments. Moreover, the rich self- 
assembling patterns including the heterochiral and homochiral phases, and their interrelationships 
are entirely governed by the same mechanism. 

PACS numbers: 81.16.Dn, 82.30.Rs, 68.43.Bc, 68.43.Hn 



Self-assembly (SA) of molecules on a solid surface is 
particularly interesting for both fundamental research 
and potential applications in the mass fabrication of 
nanoscale devices [ll-[l3|. This approach is a promising 
route to construct the functional systems with nanome- 
tre dimensions by autonomous ordering of molecules on 
atomically well-defined surfaces. The comprehensive un- 
derstandings of the mechanisms controlling the SA phe- 
nomena can steer the SA and growth processes. Vari- 
ous molecular patterns have been observed in scanning 
tunneling microscopy (STM) for organic molecules on 
metal surfaces, e.g. self-assembled supermolecular clus- 
ters and chains on Au or Ag surfaces [3|-[6|; chiral metal- 
organic complexes on Cu(OOl) Q; and a homomolec- 
ular two-dimensional honeycomb network at a Cu(lll) 
surface These organic molecules are usually pla- 

nar and rigid, and directly observable with STM. For 
biomolecules, however, due to their complicated struc- 
tures and properties such as chirality, it is not possible 
to reach an atomic resolution in STM images, which gen- 
erally appear as bright protrusions [l2|,[l3|. For example, 
recent STM studies on glycine monolayer on Cu (001) dis- 
played a blurred profile of molecules even at a tempera- 
ture as low as 5 K [l5|. Therefore, it is great desirable to 
investigate the basic molecular self-assembled structures 
beyond the STM results. 

Glycine, the simplest amino acid comprising only an 
amino group (as the head) and a carboxyl group (as 
the tail), is one of most fundamental components for 
biomolecules. Among all amino acids glycine is the only 
one that does not have chirality. However, by adsorb- 
ing on Cu surfaces glycine deprotonates and changes to 
glycinate ion. Thus, the amino and carboxyl groups 
of glycinate are strongly positively and negatively po- 
larized, respectively, upon surface adsorption, exhibit- 



ing two different chiral configurations. Plenty of struc- 
tural patterns including the p(2x4) heterochiral phase, 
c(2x4) homochiral phase, and mixed phases with differ- 
ent domain boundaries have been observed for glycine on 
Cu (001) To explain the stability of these various 

phases, Kanazawa et al. proposed a two-step mechanism, 
i.e., the c(2x4) phase is stabilized by the formation of the 
(310) steps and the p(2x4) structure by couplings of the 
molecular layer to Cu substrate [l6[ . This model ignores 
the roles of internal factors such as intermolecular hy- 
drogen bonding (HB) between glycine molecules, which 
would be the intrinsic driving force for SA. Scanning tun- 
neling spectroscopy (STS) further revealed anisotropic 
free-electron-like dispersions with electron effective mass 
differing by 10-fold along the [110] and [110] directions. 



respectively p^. Although with many tries 
building an atomistic model of glycine monolayer struc- 
ture that is thermodynamically stable and reproducing 
this anisotropic electronic behavior has not been suc- 
cessful. The mechanism governing the self-assembly of 
glycine on Cu, which has been a long-term controversy, 
is still unclear. 

In this work, we employed state-of-the-art first- 
principles calculations and extensive molecular dynam- 
ics (MD) simulations to investigate the SA structure of 
glycine on Cu(OOl). We discovered that the structure 
satisfying a complete polarity compensation reproduces 
well all known features observed in experiments includ- 
ing STM images and 10-fold anisotropic free-electron-like 
dispersion, and further predicted all possible homochiral 
and mixed phases for glycine on Cu(OOl). Our results 
indicate that molecular polarity and finite temperature 
play vital roles in determining the microscopic mecha- 
nisms for glycine SA on a metal substrate. The "ground- 
state structure" at K is not stable at room temperature; 
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instead, thermal fluctuations weaken molecule-substrate 
interaction and favor intermolecular hydrogen bonds as 
determinative forces for glycine self-assembling on Cu 
surface, which is a unique intrinsic nature of self-assembly 
in this system. 

Structure optimizations were performed within the 
framework of Density Functional Theory (DFT) using 
Vienna ab initio simulation package (VASP) [ill, [2o|, in 
which projector- augmented wave pseudopotentials and 
generalized gradient approximation [2l[ were chosen. A 
plane wave basis set was used to expand the Kohn-Sham 
orbitals with a 400 eV kinetic energy cutoff. Glycine 
molecules were placed on a seven layer Cu(OOl) slab in a 
superceh with dimensions of 10.28Ax 10.28Ax25A. The 
Monkhorst-Pack scheme [22| was adopted for the Bril- 
louin zone integration, and we have tested that a 4x4x1 
A:-point mesh is sufficient to ensure a good convergence 
in the total energy differences. The atoms in the top 
four Cu layers were allowed to fully relax until the forces 
on them were all smaller than 0.03 eV/A. Energy con- 
vergence for the geometry optimization was better than 
0.1 meV per atom. The time step in molecular dynam- 
ics simulations was chosen to be 0.5 /s. An equilibration 
process was flrst performed by slowly heating the systems 
from K to 310 K. Then we examined thermal oscilla- 
tions in atomic structure and free energy at 310 K in 
canonical ensemble using the Nose algorithm. 

Our flrst-principles calculations show that a single gly- 
cinate molecule and its enantiomeric isomer binds to the 
Cu (001) surface with an energy of 2.33 eV through one 
nitrogen atom and two oxygen atoms, all on top site of 
substrate Cu atoms, as schematically shown in Fig.dja). 
Then we examined the atomistic geometry of glycinate 
SA monolayers. Based on physical consideration and the 
characters of glycinate molecule, we found that the stable 
p(2x4) structures have the following features: First, due 
to strong attractions between the positively (amino head) 
and negatively polarized ends (carboxyl tail), glycinate 
molecules arrange themselves into a linear chain via head- 
to-tail HE along the [110] direction; Second, these chains 
form an alternating anti-parallel pattern to eliminate the 
net polarity of the two-dimensional (2D) island; Third, 
the neighboring chains are also connected by weaker HB 
along the [110] direction to make the monolayer structure 
an integral network. The above three polarity compen- 
sation reasons emphasize that the constructing glycinate 
SA arrangement is energetically favorable by maximiz- 
ing intermolecular attractive interactions while remov- 
ing the net polarity of 2D network on a surface. Under 
these restrictions, we exclusively reached a heterochiral 
p(2x4) structure of glycinate on Cu(OOl), as shown in 
Fig. [TJb). Glycinate molecules naturally alternate their 
chirality from one chain to the next, as required by the 
full polarity compensation, otherwise, less stable and iso- 
lated double-rows would form instead of the 2D mono- 
layers as observed in experiment. 
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FIG. 1: (a) Adsorbed glycinate molecule (R) and its enan- 
tiomeric isomer (^S*). Two positive sites in amino group and 
two negative sites in carboxyl group are circled by red and 
blue dashed lines, respectively, (b) The anti-parallel molec- 
ular arrangement of a heterochiral p(2x4) structure. Each 
molecule interacts with four neighbors via HB. (c) Schematic 
plot of anti-parallel HB ladder-like network. Horizontal bold 
lines represent strong HB along [110] and vertical thin lines 
weak HB along [110]. The arrows indicate the HB directions. 



Geometry optimization shows that all 0...H lengths in 
the HB along [11] are 1.6 A with the binding energy of 
210 meV (indicating a strong HB), whereas those along 
the [110] direction are 2.0 A with the energy of 138 meV 
(weak HB) [2^. They are denoted in blue and green 
dashed lines in Fig. [TJb), respectively. This anti-parallel 
HB network is also schematically drawn in Fig. [U^c), to 
indicate the arrangement of molecular polarity. In this 
structure each glycinate acts as two HB donors and two 
HB acceptors^ where all its polar sites are nicely satu- 
rated with HB [24]. This is consistent with the general 
understanding of H bonding interactions: namely, the 
selectively, linearity, and saturation of H bonds. 

Any structures proposed for glycine SA on Cu(OOl) 
would have to withstand stringent tests to reproduce 
properties measured experimentally; our structure model 
based on polarity compensation, shown in Fig. [TJ satis- 
fles these tests very well. To conflrm it is indeed the 
one observed in experiment, we compare STM images, 
electronic band structure, structural phases, and ener- 
getics from both theory and experiment. Fig.[2fa) shows 
the simulated STM image for the p(2x4) phase, together 
with the measured image [l5| in the inset for compari- 
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son. The STM image is calculated by integrating lo- 
cal states with energy <100 meV below the Fermi level 
within Tersoff-Hamann approximation [25[. It displays a 
pattern of triangular protrusions: triangles slightly tilted 
to the left are glycinate in S conformation, while those 
tilted to the right are in R conformation, forming an 
array of alternating anti-parallel rows. In experiment a 
pattern of blurred triangles with two orientations was ob- 
served, forming anti-parallel alternating rows along [110], 
in a very similar manner to simulated images. Within 
each row the molecules are close to each other, indicat- 
ing a stronger interaction within the rows; while there is 
a larger separation between the rows, suggesting weaker 
row- to-row interactions. These features are well repro- 
duced and explained by our model. 

More stringent test comes from the 10-fold difference 
in electron effective mass m* for the free-electron-like dis- 
persions along the [110] and [110] directions observed 
in STS Again this strong anisotropic behavior is 

obtained for our proposed p(2x4) structure. We calcu- 
late the surface band structure of glycinate on a single 
Cu(OOl) layer, in order to remove the noisy background 
of bulk substrate. This approximation does not affect 
the anisotropic property because the removed Cu atoms 
are symmetric to the [110] and [110] directions. The re- 
sults are shown in Fig. [2fb). There is only one band 
across Fermi surface for both [110] and [110] directions. 
The solid lines represent the calculated results, and the 
dashed lines are rescaled by multiplying a constant 2.46 
for direct comparison with experimental data. Dashed 
curves perfectly match the experimental results, show- 
ing parabolic free-electron-like dispersions, as well as the 
10-fold difference in the effective mass with respect to 
the [110] and [110] directions. The scaling constant is 
attributed to an enlarged effective mass m* by using a 
single layer Cu in our approximation. Meanwhile, the 
Fermi energy is also shifted to include contributions from 
Cu layers underneath. Clearly, the anisotropy originates 
from different strengths of interactions between neighbor- 
ing molecules in two directions, i.e., strong HB along the 
direction [110] and weak HB along [110], shown schemat- 
ically in Fig. ffl^c). Note that the weaker HBs are mis- 
aligned along [110] which would significantly decreases 
the conductivity in this direction as well. Thus, different 
HB strengths together with their arrangements result in 
the 10-fold difference in the effective mass. 

To interpret their STM observations, Kanazawa et al. 
adopted a structural model consisting of only head-to- 
tail alignments in the same direction, which is prevailing 
in literature for glycine SA on Cu(OOl) [13, [ll, E [26| . 
For convenience of presentation, we name this struc- 
ture the head-to-tail (HT) model. This model, how- 
ever, does not consider the influence of polarity to side 
neighbors, not mention to the HB network. Thus, nei- 
ther STM images nor free-electron-like dispersion in ex- 
periments were reproduced by the HT model. Indeed, 
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FIG. 2: (a) The calculated STM image for the p(2x4) struc- 
ture, the insert is experimental result from Ref. [15] for com- 
parison, (b) Energy dispersion relations for the [110] direction 
in blue and [110] in red, respectively. Solid lines represent our 
calculated results and dashed lines are rescaled for comparison 
with experimental results from Ref. [15] (blue dots for [110] 
and red triangles for [110] ). The Fermi energy is indicated 
by horizontal dashed line. 



first-principles calculation shows that the HT structure is 
76 meV higher per molecule in energy than our structure 
for free-standing monolayers, indicating that our model is 
intrinsically more stable for SA than the HT model. In 
our p(2x4) structure glycinate molecules are saturated 
with H bonds and all H bonds conform the ideality of 
H bonding interactions such as selectively and linearity, 
thus its free-standing structure is more stable than the 
unsaturated HT structure. The apparent -CH...0 hydro- 
gen bonds with the length of 2.3 A in the HT structure 
cannot be regarded as an ideal HB. 

However, a complexity arises from the subtle bal- 
ance between the intralayer molecule-molecule interac- 
tions and the molecule-substrate binding, which is a cen- 
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FIG. 3: Free energy fluctuations as a function of time in 
molecular dynamics simulations at 310 K for our structure 
and that adopted by Kanazawa et al. (Refs. [15,16]). The 
inset shows the distribution of free energies for the two cases. 



tral theme in all SA processes. At K, first principles 
calculations show that the HT model is lower in total en- 
ergy by 140 meV per glycinate than the adsorption of our 
structure based on polarity compensation. That is, the 
HT model is the "ground state structure" for glycinate 
layer on Cu(OOl) at K. Then why it is not observed in 
experiments, which shows very different STM/STS fea- 
tures as to those of the HT structure? 

To solve this controversy, we performed ab initio MD 
simulations for both HT and our structure model at room 
temperature (RT)~310 K. Both configurations are stable 
in trajectories lasting for 5 ps. However, we found the 
"ground state structure" at K presents a higher free en- 
ergy at 310 K, which is now 75 meV per molecule higher 
in energy than our model, shown in Fig.[3l The difference 
in free energy is coincidental to that for free-standing lay- 
ers, 76 meV, suggesting that the SA structure at RT is 
dominated by intermolecular H-bonds, rather than the 
molecule-substrate bonds as in the case at K. At abso- 
lute zero the unsaturated HT model has stronger interac- 
tion with the substrate compared to our saturated struc- 
ture. Thermal fluctuations weaken glycinate-Cu bonds, 
and free glycinates adjust themselves to maximize HB 
interactions while maintaining glycine- Cu interaction in 
a similar strength. 

In fact, besides the molecule- molecule and molecule- 
surface interactions, the temperature plays as the third 
important parameter in molecular self-assembly at sur- 
faces During the SA process, molecules interact with 
each other sufficiently with the help of thermal energy. In 
experimental conditions (temperature ~ 370 K[15| and 
~ 430 K[12]), glycinate molecules move around on Cu 
surface in translational and rotational modes, and inter- 
act with each other by HB. The HB would break and 



re-bond, providing sufficient opportunities for molecules 
to search for stronger intermolecular interactions. This 
process, together with strong polarity compensation of 
glycinate, promotes formation of anti-parallel HB net- 
work on Cu(OOl). Once formed, the structure is trapped 
as a metastable phase even at very low temperatures. 
The fact that the structure we proposed is observed at a 
temperature as low as 5 K in experiments, implies that 
transition from our structure to the "ground-state" HT 
structure could be blocked by a large energy barrier. We 
emphasize that the SA structure at a finite temperature 
is determined by the free energy, rather than the total 
energy at K. The molecule-surface interaction decides 
the adsorbed glycinate molecular configuration in the tri- 
dentate fashion, but the intermolecular interaction plays 
a dominant role for the glycinate SA arrangements on 
Cu(OOl). Therefore, it is demonstrated that at finite 
temperature the SA on surfaces are mainly determined 
by the relatively weak HB rather than the stronger in- 
terfacial interaction. 

Based on above investigations, we summarize a gen- 
eral rule for the construction of glycinate SA structures 
on Cu surface. The rule has three elementary compo- 
nents: (i) Because of the strong polarity, the positive site 
of a glycinate interacts with the negative site of another 
via hydrogen bond to form a head-to-tail chain, (ii) The 
polarity of glycinate also affects its side neighbors, i.e., 
the positive sites of neighbors will approach to the neg- 
ative site of the glycinate, and vise versa. Thus, (i) and 
(ii) together result in the local molecular configurations 
where the positive site of each glycinate is surrounded by 
the three neighbors' negative sites, and the negative site 
surrounded by the three positive sites, (iii) Each glyci- 
nate molecule contributes its all active sites, two positive 
sites in amino group and two negative sites in carboxyl 
group, to interact with its neighbors via HB. According 
to the rule (iii), glycinate molecule will automatically 
adjust its chirality, an important feature associated with 
SA, to self-organize HB network structures. These three 
principles for constructing glycinate SA are collectively 
called the polarity compensation rule here. 

Besides the dominant heterochiral ^o(2x4) phases 

homochiral 



12 



widely observed in experiments 
phases were also found by several groups for glycinate 
SA on Cu(OOl) The presence of both heterochi- 

ral and homochiral phases induces a variety of SA pat- 
terns in STM images [l6|, e.g. a homochiral array at the 
boundary between heterochiral p(2x4) phases. To reach 
a unified understanding on these rich SA patterns, we 
invoke the polarity compensation rule to unravel the un- 
derlying mechanism for the formation of different chiral 
phases. 

We first look at the boundary between two p(2x4) 
domains, pi(2x4) and p2(2x4), where p2 is shifted by 
one lattice of Cu surface along [110] with respect to pi, 
shown in Fig. HI At the boundary, the glycinate with 
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FIG. 4: Molecular arrangements for comparing with STM 
images in Ref. [16], where exhibit various structures including 
the heterochiral j9(2x4) and homochiral 0^(2x4) domains and 
their relationship. Dashed lines indicate domain boundaries. 

R chirality has to turn its head into S chirality in order 
to form a strong HB across the boundary and a weak 
HB with the side neighbor. Thus all molecules at the 
leftmost column of p2(2x4) become S chirality, denoted 
as pg{2x4:) in Fig. HJ This homochiral glycinates array 
must exist between the shifted heterochiral pi and p2 
domains, offering full HBs to stabilize the boundary and 
lowering the energy by ~100 meV per molecule. This 
simple mechanism provides a natural explanation for the 
well observed homochiral phase between p(2x4) domains 
in experiment, indicated by white arrows in Fig. 2(e) of 
Ref. [16]. 

The above observation demonstrates how glycinate 
changes its chirality to interact with its neighbors ac- 
cording to the rule when the molecule shifts along [110] 
by a Cu position. Shifting continuously in this way, all 
glycinates with R chirality would change to S chirality in 
order to bond together, so a complete homochiral phase 
of S chirality is formed. More importantly, this homochi- 
ral phase naturally develops a step along the [310] direc- 
tion, on which all molecules are saturated with HB and 
the step becomes inactive, shown in the left corner of 
Fig. HJand Fig. [5l This triangular p|^(2x4) phase is in- 
deed observed by STM, denoted as c(2x4) in Fig. 2(e) 
of Ref. [16]. 

Since glycinate is enantiomeric, it is expected that an 
alternative homochiral c^(2x4) phase with all R chiral- 
ity also exists. Indeed this homochiral i?-phase is formed 
in the same way at the right side of the p(2x4) pat- 
tern, illustrated in Fig. [5l Similarly, this phase could 
also be formed by shifting the glycinate molecules along 
[110] at the left side of the p(2x4) pattern by a Cu posi- 
tion, while the c|^(2x4) phase at the right side. Both c'^ 
and c'^ phases are composed of hydrogen-bonded twin 
chains. Besides the aforementioned homochiral phases, 
by shifting the glycinate a Cu position along [110] se- 
quently, other two homochiral phases above and below 



the p(2x4) pattern can also be constructed based on our 
rule, named as C5'(2x4) and Ci^(2x4) in the upper and 
lower part of Fig.O respectively. The chirality of the out- 
ermost glycinate in the p(2x4) phase decides the chirality 
of the c phase. Since each adsorbed glycinate molecule 
occupies a (2x2) Cu surface cell, every possibility for the 
c' and c phases by shifting a Cu position along each di- 
rection has been considered. It should be noted that all 
above homochiral phases are different from the c(2x4) 
phase model adopted by Kanazawa et al.jld]. 

We have exhausted all possible SA structures of gly- 
cinate on Cu(OOl): a unique heterochiral p(2x4) phase 
and four homochiral phases based on our rule. In par- 
ticular, the chirality of each phase is automatically de- 
termined by the rule (iii). The interrelationships be- 
tween various phases are displayed in Fig. [5] where the 
heterochiral p(2x4) phase is surrounded by different ho- 
mochiral phase in each direction. The reason that we 
put them together is to emphasize their seamless con- 
nections where the boundaries between the hetero- and 
homo-chiral phases also obey our rule. No matter where 
a glycinate molecule locates in single phase or bound- 
ary, it fully interacts with the neighbors by two positive 
sites in amino group and two negative sites in carboxyl 
group via HB without any exception. To our knowledge, 
glycine on Cu(OOl) exhibits the richest 2D SA structures 
on surface, in which all SA patterns and their interre- 
lationships are naturally explained under a single mech- 
anism. Since p(2x4) and c|^(2x4) structures have al- 
ready been observed in STM experiment [l^ , searching 
for other predicted phases experimentally remains a chal- 
lenge to prove our theory. 

To conclude, we have clarified that the self-assembly 
structure is determined by the free energy at finite tem- 
perature, and revealed the mechanism for self-assembly of 
glycine on Cu(OOl) based on full polarity compensation. 
The ideality of H bonding environment reaches a unique 
antiparallel head-to-tail network of glycinates in p(2x4) 
phase. We elucidate that thermal fluctuations at finite 
temperatures weaken the molecule-substrate interaction 
and favor the intermolecular interactions, which plays 
a dominant role to orchestrate glycine self-assembly. 
The long-term controversy about competing phases, i.e., 
homechiral phase and heterochiral phase, as well as their 
relationships has been uniformly resolved in the same 
framework. The present work represents a distinct ex- 
ample that self assembly of molecules is driven by in- 
termolecular interactions themselves, rather than being 
fixed by interacting with templates, even though the 
molecule-template interaction is stronger at K. Al- 
though demonstrated for a small molecule on a simple 
surface, the SA mechanism and the rule revealed here 
are expected to have much wider implications in under- 
standing self-assembly of polar molecules with other tem- 
plates and of larger biomolecules in biologically relevant 
environments [1,1H. 
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FIG. 5: Panorama of the SA monolayer of glycinate on Cu (001). There are four homochiral phases around the heterochiral 
phase j9(2x4). Domain boundaries are represented by the dashed hues, and the unit ceh for each phase is denoted by sohd 
hues. For the sake of clarity, the Cu substrate is not shown. 
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